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1.  Introduction. 


Recently,  N.  N.  Bojarski  (1974)  conjectured  that  inverse 
source  problems  for  the  wave  equation  could  be  investigated 
by  means  of  a Fredholm  integral  equation  of  first  kind  which 
took  the  form, 


(1.1) 


/dK(1  - £*)  P(l*)  d^r*  = 0(r), 


Here,  p(r^')  denoted  the  unknown  source  whose  support,  however, 
was  a^jumed  a priori  to  be  located  within  a domain  D.  Further- 
more, 0(0  denoted  a function  determined  by  observations  of  the 
field  on  the  boundary  of  D. 

Later  that  year,  N.  Bleistein  and  N.  N.  Bojarski  (1974) 
presented  an  invited  talk  at  the  Summer  Institute  on  Inverse 
Problems  at  the  University  of  California  at  Irvine,  in  which 
they  gave  a rigorous  derivation  of  (1.1),  identifying  the 
kernel  K as  the  difference  of  the  free  space  outgoing  and 
incoming  Green's  functions.  The  equation  as  thus  formulated 
has  a unique  solution,  but  it  was  soon  realized  that  the  source 
function  p,  which  arose  from  a time  transform  of  the  wave 
equation,  would  in  general  depend  on  the  transform  variable, 

0).  Thus,  the  treatment  and  numerical  results  presented  by 
these  authors  only  applied  to  the  case  in  which  the  time- 
dependent  source  had  the  form,  p(£)  6(t),  where  6(t)  represents 
the  Dirac  delta  function. 

Accordingly,  the  present  authors  (1975)  gave  a derivation 
of  the  integral  equation  valid  for  general  time-dependent  sources 
and  showed  that  in  this  general  case,  the  source  cannot  be  deter- 
mined uniquely.  We  gave  several  characterizations  of  the 


» 


m 


p 


non-uniqueness  and  also  gave  illustrations  of  some  type  of 


additional  a priori  information  about  the  source,  which  would 


allow  its  unique  determination.  Furthermore,  we  extended  the 


theory  to  the  Maxwell  Equation  system. 


In  a sequel,  now  in  preparation,  we  will  present  additional 


examples  of  both  the  source  determination  and  source  synthesis. 


Many  of  these  applications  proceed  most  easily  by  treating  the 


spatial  Fourier  transform  of  the  integral  equation.  However, 


we  decided  that,  since  the  integral  equation  is  of  quite  clas- 


sical type,  a study  of  its  eigenfunction  structure  from  the 


Hilbert  Space  viewpoint  would  be  illuminating  and  it  is  that 


study  we  present  here. 


2.  Notation  and  Derivation  of  the  Integral  Equation. 

We  deal  with  functions  of  space  and  time  denoted  by 


F(r,t),  where 


(2.1) 


£=(x,y,z),  r = I r i , r = r/r. 


Tne  time  transforms  of  such  functions  are  denoted  by  the 
corresponding  lower  case  letters. 


(2.2) 


+“> 

f(£,aj)  = / e^“^F(£, t)dt 


and  we  shall  occasionally  also  have  reference  to  the  time- 
space  transform,  denoted  by  a tilde  as, 


(2.3) 


f(k,0))  = ///f(£,o))e  ^^r. 


We  consider  the  wave  equation, 


(2.4) 


(V2  _ ) u(r,t)  = -F(r,t) 


and  assume  that  we  know  a priori  that  the  source  F has  spatial 
support  within  a sphere  of  radius  a and  also  that  the  source 


has  temporal  support  in  the  interval*,  ^-t  < t < t . We  assiune  that 

0 0 


(2.5) 


U(r,t)  E 0 for  t < -t 


Upon  time  transformation,  we  obtain 


(2.6) 


[V^  + (^)2J  u(r,u))  = -f(r,u) 


with  u satisfying  the  outgoing  radiation  condition. 


(2.7) 


1 — i— r '' 

u(r,m)  - — e c u (r,o))  as  r 

4rr  0 


* This  can  be  weakened  somewhat. 


V, 


■■ — 


Here,  u is  a quantity  analogous  to  the  scattering  cross  section 
0 

called  the  phase  and  range  normalized  scattering  amplitude. 

We  now  present  a simplified  derivation  of  the  generalized 
Bojarski  integral  equation  which  proceeds  immediately  from  the 
Green's  Identity, 


(2.8)  ///p  |u(£',(o)  Lv(^',^,o))  - v{^' ,2^,0.')  Lu(r_',u)|  d^r' 


“ //  fuV'v  - vV'ul  ds' 

3D 


(2.9) 


L = V’2  + (i!i)2, 
c 


3D  denotes  an  (observation)  boundary  which  lies  outside 
the  sphere  r £ a (the  support  of  the  source) . D is  the  interior 

A> 

of  3D  and  n'  is  the  unit  outward  normal  to  D.  We  choose  v as  a 
solution  of  the  homogeneous  Helmholtz  equation. 


(2.10) 


[v2  + (§)2]  v(r,r',o))  = 0 


which  is  regular  at  its  source  point  £ = £* . The  simplest 
choice  of  v is 


(2.J1) 


V = j/f  R),  R = |r  - r’l; 


ioW  = 


sin  X 


Here,  j.  is  the  spherical  Bessel  Function  of  first  kind  and 


order  zero. 


On  using  (2.10)  and  (2.6)  in  the  Green's  Identity  (2.8), 


we  have  at  once  the  integral  equation 


(2.12) 


where 


(2.13) 


(2.14) 


Kf  = 0(£,O)) 


Kf  = ///J  (— ) f(r'.o))  d^r' 
■’•’•’D  0 c — 


0=//  r>  • [u(r',co)Vj  (“I) 
•’■'3D  0 c 


■ 3„(^)  v'u(_r’ ,(0)]ds' . 

0 (- 


The  integral  equation  (2.12)  is  valid  for  all  £,  and  since  f 
has  support  in  a sphere  of  radius  a,  the  integration  in  (2.13) 
can  be  regarded  as  extending  either  over  all  space  or  just  over 
the  a-sphere  as  we  find  conceptually  convenient. 

The  integral  equation  (2.12)  differs  only  slightly  from 
Bojarski's  original  conception.  The  kernel  and  0 differ  only 
by  a common  multiplicative  factor  from  his,  and  most  importantly 
0 is  still  an  observable  quantity  computed  from  observations  of 
u and  3u/9n  on  the  single  surface  3D. 


As  mentioned  in  the  introduction,  most  of  the  practical 
applications  of  the  integral  equation  follov?  most  easily  from 
its  spatial  Fourier  Transform,  which  we  remark  in  passing,  is 


given  by 


(2.15) 


1c(k,ck)  = 0(k) 


where 


(2.16)  0(k)  = -//  e“^—  - n'  • [iku(r',ck) 


+ V’u(_r*  ,ck)1ds’  'v-  u (k,ck)  as  r 

0 


3.  Eigenfunction  Analysis  of  the  Integral  Equation. 
We  study  the  problem 


(3.r< 


= Xij;(_r,w),  |£|  ^ a,  |_r ' | £ a, 


where  K is  defined  by  (2.13).  and  is  a ‘'‘P-gutar  function. 

Since  the  kernel,  - _r ' | ) , of  K is  a symmetric,  continuous 

function  on  a compact  product  space,  we  are  dealing  with  the 
most  classic  case  of  symmetric  Hilbert-Schmidt  theory.  Below, 
in  statements  1 /,  we  shall  summarize  the  principal  conclusions 
o'  jis  theor/  and,  when  appropriate,  give  the  specializations 
- r these  results  to  our  kernel.  Note  that  our  eigenfunctions 
exist  only  in  the  a-'--phere  and  that  the  precise  value  of  a is 
somewhat  arbitrary.  Also,  we  assume  that 


(3.2) 


ip  (r,o))  c - [a]  , 
— 2*- 


where  the  notation  indicates  that  ip  is  square-integrable  on  the 
a-sphere. 

We  begin  our  review  of  the  theory  with  the  following 
definitions . 

1.  E = linear  manifold  generated  by  the  eigenfunctions 
for  X 0. 

N = null  space  of  K (i.e.,  the  solutions  of  Kf  = 0) . 

3.  E and  N are  subspaces;  in  fact,  L [a]]  = E ® N.  This 

result  means  that  for  any  f e L [al , 

2 

4.  f = Ef  i{)  + h , where  the  f are  the  unique  Fourier 

n n n ^ 

Coefficients  of  the  eigenfunctions  which  generate  E, 
while  h is  an  arbitrary  element  of  N. 


£"J2tS^S?SSS'a2steo« 


For  our  kernel,  the  eigenfunctions  generating  E are 


explicitly  given  by 


(3.3) 


"^0  (6.4>). 

«,m  — No  *.  c £m 


£ = 0,1,2  ...  ; |m  < £ 


The  corresponding  eigenvalues  are  given  by 


(3.4) 


A = A = 4tiN  2 
£m  £ £ 


and  the  normalization  factor,  N^,  is  given  by 


(3.5) 


N,2(<^,a)  = j,2(“I)r2dr 


To  establish  (3.3),  recall  that  the  kernel,  j , was  chosen 

0 


as  a regular  solution  to  the  homogeneous  Helmholtz  Equation.  Thus, 


upon  applying  the  operation,  + (^)^,  to  both  sides  of  (3.1), 


we  have  at  once 


(3.6) 


a[v2  + (^)2]^  = 0 


In  determining  E,  we  assume  A =j=  0,  so  it  is  clear  that 


eigenfunction  in  E is  a regular  solution  to  the  homogeneous 


Helmholtz  equation.  If  we  expand  ij>  in  the  complete  set. 


Y„  (6,4>)»  of  regular  spherical  harmonics  as 
Jom 


oo  z 


(3.7) 


i|)(r,o))  = Z E f (r,u)Y  (6,4)), 
£=0  m=-£ 


we  find  that  f^^(r,u)  satisfies  the  second  order  spherical 


Bessel  or\  'nary  differential  equation.  The  only  L [a^ 


v*'}-  -•■  c 

s, 

'ShV'- 


«tt^‘  ••  '-»!  v’ 

I® 

t#-- 

fftri,h-'  - ' 


iff-' 

wM: 
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&> 
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solution  is  j (^  r) , so  that  « j ^ 


concludes  the  proof  of  (3.3), 


If  we  expand  our  kernel  in  spherical  harmonics  as 


00  I 


(3.8)  ,i  (^  ir  - r’!)  = 4x  Z Z j (iill) j (“ll) 


i!.=  0 m=-£ 


£ c Z c 


4>„  (e.4>)t>„  * (0,({>), 

ilm  «-m 


(see  Jackson,  1962),  it  is  easy  to  compute  Ki|)  , and  hence 

xm 


obtain  the  eigenvalues  A„  given  in  (3.4). 

x.m 


To  any  function  f in  L^[aJ,  we  can  associate  both 


an  expansion  in  spherical  harmonics  and  an  eigen- expansion. 


To  distinguish  these,  we  henceforth  use  f^^^^  to  denote  the 


spherical  harmonic  coefficients  and  f„  to  denote  the  eigen- 

x,m 


coefficients.  Thus, 


00  n 


(3.9) 


f(r,o))  = Z Z f„  (r,(o)Y„„(e,^) 


fc=0  m=-it 


(3.10) 


00  Z 
x.=  0 m=-2. 


The  relation  between  these  coefficients  is 


(3.11) 


(3. .2) 


f.  (r,(o)  = //  f(r,u))Y„  *(0,<^)dn 


(3.13) 


(“)  = ///  f(r,u)'|'  *(r,o))r2drdJl. 
X,m  “ Jem  ~ 


Here,  dQ  denotes  an  integration  over  the  unit  sphere. 


The  general  expansion  result  (4)  can  now  be  explicitly 
expressed  for  our  kernel  as 


(3.14) 


f ^ ^ hmhm  + 

Z=0  m=-Z  ^ 


where  the  f are  the  unique  eigen-coefficients  (3.13),  and  h 
Jim 

is  an  arbitrary  element  of  N. 

We  now  turn  to  the  characterization  of  the  range  of  K. 


We  write 


(3.15) 


and  have 


H(K)  = range  of  K, 


5.  T1(K)  is  a linear  manifold  in  E,  which  is  not  closed. 
Furthermore,  we  have 

6.  geH(K)  ->  the  eigen-series  associated  with  g is  uniformly 
convergent  (instead  of  merely  convergent). 

This  latter  result  requires  that  the  iterated  kernel  be 
bounded  and  continuous,  which  is  certainly  true  here  since  we 
deal  with  a continuous  kernel  on  a compact  set.  The  second 
part  of  (5)  requires  only  that  the  kernel  be  Hilbert-Schmidt 
and  there  exist  infinitely  many  eigenfunctions  (i.e.,  that  K 
be  non-degenerate) . 

The  fact  that  *R(K)  is  not  closed  means  that  the  alternative 


theorem  does  not  apply,  but  since  N = always  holds,  we  do 


of 


E = closure  (R) . 

Using  our  eigenfunctions,  we  see  at  once  that  the  solutions 


11 


(3.16) 


Kf  = g(^,i>)),  geR(K), 


r < a 


are 


(3.17) 


f(r,o))  = I 


I i'H*  . , 


£m 


a=o  m=-ii  i. 

where  h is  an  arbitrary  element  of  N.  Thus,  the  range  may 
be  characterized  by 

£ gj 

^ r 

2.=  0 m=-ll  '"Z 


(3.18) 
or  by 

(3.19) 


geR(K) 


E E ^ eL  [a] 
X /.  )lm  2 ^ 


ge11(K) 


00  Z 
E E 
Z=o  m=-£ 


’i.m 


< + o> 


where  the  g are  the  eigen-coefficients  of  g.  The  existence 

2.m 


criterion  (3.19)  is  quite  stringent  for,  while  it  is  also  true 
for  non-degenerate  symmetric  compact  operations  that  the  sequence 


of  eigenvalues,  X , satisfy  X -»•  0 as  n -»■  “>  , we  have  for 

n n 


our  operation  extremely  rapid  convergence  to  zero.  In  fact,  the 
explicit  expressions  (3.4)  and  (3.5)  for  the  eigenvalues  yield 


(3.20)  = 5^(|t|)'’'Ll  + 0(|)]  . » 


Thus,  in  a sense,  ^(K)  is  much  smaller  than  E. 

In  the  following  table,  we  give  some  typical  values  for 


X /2ira^  so  that  the  reader  can  judge  the  noise  amplification 


inherent  in  an  eigenfunction  solution  such  as  (5.17). 


i 


12 


\ coa 
«,  \ c 

.1 

1 

10 

100 

0 

.6653 

.5454 

.9544x10-2 

— 

1 

.4438x10~3 

.3850x10-1 

.1040x10-1 

— 

2 

.A122xl0~‘^ 

.1136x10-2 

.9174x10-2 

— 

3 

.2014x10“^° 

. 1 840x1 0"*" 

.9452x10-2 

— 

5 

.1423x10-16 

.1332x10-6 

.7783x10-2 

. 9950x10" 

10 

— 

— 

.6130x10-2 

.9943x10"'' 

15 

— 

— 

.8437x10-2 

.9929x10"'' 

Table 

of  values  for 

* _ ^ 2/Oia'. 

~ It. 

27ra6 

- -i 

(^) 

W=  conclude  this  section  by  providing  several  characterizations 
of  the  null  space  N . Aside  from  the  obvious  characterizations, 

Kf  = 0 and  feE-*-  , we  have 


(3.21) 

(3.22) 


= 0;  J.  = 0,1,2,  ...  ; m £ X,, 


feN  f(k,ck)  = 0, 


(3.23) 


feN  f is  a non-radiating  source 


By  the  phrase  "non-radiating  source",  we  mean  one  whose  field, 
outside  the  source  region  r _<  a , vanishes  identically.  This 
last  characterization  identifies  the  mathematical  concept  of  a 
null  space  with  those  sources  which  are  physically  undetectable 
outside  the  source  region.  The  proofs  of  statements  (3.21)-(3.23) 
appear  in  Bleistein  and  Cohen  (1975).  Here,  we  show  in  a con- 
structive way  that  (3.21)  implies  that  N is  infinite  dimensional. 
Consider  the  functions 

( /o  f(r’,a))jj^,(^’)r-2dr*-j 

(3.24)  h(r,m)  = -^f(r,a)) ^ 2 ( Y , ,(6,<j.) 

•'o  c 


For  these  functions,  we  have 
(3.25)  hp„  = /!  / h(. 


/g  / h(_r,o))ij)^jjj(r,oi)r2dfidr 


^ Jo 


0 -"V  C 


j (— )r^dr  -tS  , 

-’ll'  c £,S.  m,m' 


U-  „ ,2ar  - ^0 

N \ 0 ^r^dr 

jC 


j ^r^dr 
% Z 


m 

_ f . 

i 

(‘il  )r'2dr’  1 

i 

c ( 

i 

dr'  / 

0 


Hence,  for  each  f(r,a))  and  each  harmonic  Jl,m,  there  is  a 
solution  in  the  null  space.  For  example,  if  we  choose  the  (0,0) 
harmonic  and  f = r*^,  k = 0,3,2,  ...  , we  obtain  an  infinite 
family  for  which  the  quadratures  can  be  explicitly  performed.  We 
also  note  that  if  f(r,u)  is  an  entire  function  of  order  one 
in  0),  so  is  h(r,(<)).  Thus,  the  inverse  transform  H(r^, t) 
would  have  compact  support  in  time  and  would,  thus,  be  a legitimate 
source  if  only  f(r,o))  vanishes  for  r > a. 


mmH 


4.  Existence  and  Uniqueness  for  the  Inverse  Source  Integral  Equation. 


We  can  now  state  the  fundamental  existence  theorem  for  the 


integral  equation  (2.12) 


Theorem  1:  The  equation  Kf  = 6 has  solutions  if  and  only  if 


(4.1)  E E I//  n'  • *V’u]ds'12  < 

j=0  m=-£'8D  ^ 

Moreover,  if  the  equation  has  solutions,  then  it  has  the  infinite 
family  of  solutions. 


(4.2)  f = h + E E //  n’  • [uV'tl)  * ~ *V’u]ds'i|)  , heN. 

£=0  m=-«,  3D 

Proof;  Kf  = 0 has  solutions  if  and  only  if  OtR(K),  hence,  by 
(3.18),  if  and  only  if 


(4.3) 


s z -jm  < CO  . 
£=0  m=-£ 


If  it  does  have  solutions,  then  by  (3.17)  they  are  given  by 


(4.4) 


« £ On 

f = h + E E ^ , heN. 

n n n A „ £m 

£=0  hf=-£  £ 


Thus,  to  establish  the  theorem,  we  need  only  show  that 


(4.5) 


Q = X //  n'»  juV’t);  * - ijj  *V'\T]ds'. 


£m  £ 


£m  £m 


This  can  be  done  by  recalling  the  definition  of  0 in  (2.14) 

and  realizing  that  the  expansion  (3.8)  for  j (—  R)  can  be 

0 c 

written  in  terms  of  the  eigenfunction  as. 


(4.6) 


0 jt=0  m=-£  ^ ^ ^ 


Using  this,  we  obtain 





00  Z 


(4.7) 


e(r,(D)  = E E //  "’  * t *(r',u)) 

— * •'•^  ■"  Jim  ~ 


Jl=0  m=-£  ^ 3D 


which  immediately  gives  (4.5)  and  establishes  the  theorem. 

We  remark  that  if  0 is  computed  from  data  from  an  actual 

source  feL  ra"]  , then  the  existence  criterion  in  Theorem  1 
2 

is  always  satisfied  because  the  integral  equation  itself 
implies  that 


(4.8) 


0 = 

£m  ^ Ji^m 


and  hence 


(4.9) 


00  5, 

E E 
)l=0  m=-Jl 


00 

E E 
£=0  m=-Jl 


< CO 


» 


since  feL  fa] . 

We  further  point  out  that  because  of  the  uniform  convergence 
of  the  eigen-expansion  of  functions  in  the  range  and  the  final 
remark  made  at  the  end  of  the  last  section,  infinitely  many 
solutions  of  the  form  (4.2)  can  be  found  which  will  give  rise 
to  physical  space-time  sources  F(_r,t)  if  the  hypothesis  of 


the  theorem  is  met. 


5.  An  Application. 

As  mentioned  in  the  introduction,  we  plan  a sequel  in  which 
various  applications  of  equation  (2.12)  will  be  presented.  Natu- 
rally, the  key  to  these  applications  is  the  restriction  of  f to 

subsets  M of  L Fal , such  that  the  problem 
2 “ 


(5.1) 


Kf  = 0 feM 


has  a unique  solution.  We  have  already  given  some  examples  of 
suitable  choices  of  M in  Bleistein  and  Cohen,  1975.  Here,  we 
confine  ourselves  to  mentioning  just  one  such  application.  Namely, 
we  point  out  that  if  0 satisfies  the  existence  criterion  and, 
furthermore,  if  we  impose  a priori  that  fcE,  then  clearly  the 


unique  solution  is 


(5.2) 


“ 0, 


f(r,o))  =11 


Z=:  m=~Z  Z 


The  restriction  fcE  is  not  devoid  of  interest,  for  by 
the  characterization  (3.23)  of  the  null  space,  any  component 
of  feE-^  does  not  radiate.  Thus,  any  energy  used  in  producing 


such  a component  is  "wasted".  In  this  sense,  the  solution  (5.2) 
is  the  most  economical  of  the  family  of  solutions  (4.2)  and,  thus, 
is  of  particular  interest  in  the  synthesis  problem. 
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